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1. Introduction 

Consider the one-dimensional gas of spin 1/2 fermions interacting via the 5-potential. 
The Hamiltonian of this model is 

H = f A J_ E , imMx) + £,„,(.)„, M - ,„( J . (1.1) 

The fermionic fields ip a (a is a spin index, a =f, j) satisfy canonical equal-time anti- 
commutational relations 

*Pa(x)ipl(y) + ipl(y)ip a (x) = 5 al3 5(x - y). (1.2) 

The operators n a (x) = 4'l l (x)4' a (x) are the density operators for spin-up (spin-down) 
fermions and n(x) = n^(x) + n^(x) is the total fermion density operator. In this paper 
we will consider the infinite U limit of one-particle correlation functions 

G h (x,t) = (^(x,t)^(0,0)) (1.3) 

and 

G e (x J t) = (^(x,tW(0 > 0)> (1.4) 
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where the average () is taken over the thermal ensemble at a given temperature T and 
chemical potential \x. We will calculate the large x and t asymptotics of the correlation 
functions (jl.3j) and (jl. 4)1 assuming that (a) the limit of infinite repulsion U = oo is 
taken at finite positive chemical potential //, and (b), the zero temperature limit is 
taken afterwards. The ground state of the system is infinitely degenerate at infinite 
U. This degeneracy leads to a non-trivial (non-conformal) behaviour of the correlation 
functions (jl.3j) and (jl.4|) at zero temperature. The physics of this phenomenon were 
presented earlier [TJ, here we present a detailed derivation of these results. 

Our starting point is the determinant representation for Gh(x,t) and G e (x,t), 
recently derived by Izergin and Pronko [2|. We relate the objects entering this 
representation to the solution of the corresponding matrix Riemann-Hilbert problem 
(RH problem) jHj . We perform the asymptotic analysis of this RH problem making use 
of previous work 01 Ej. The large x asymptotics of the correlation functions at t = are 
given in subsection 14.31 The asymptotics of the time-dependent correlation functions 
in the space-like region are given in subsection 16.31 The asymptotics in the time-like 
region for x/t = const > are given in subsection 18.21 The large t asymptotics of the 
correlation functions for x = are given in section El 

2. Determinant representation and the Riemann-Hilbert problem for the 
correlation functions of the model 

This section is introductory. In subsection 12.11 we recall some elementary relations for 
the correlation functions (jl.3j) and (jl.4j) . In subsection 12.21 we recall the results of 
[2] on the determinant representation for these correlation functions. In subsection 
12.31 we examine the analytic properties of the functions entering this determinant 
representation. In subsection 12.5) following we formulate the matrix RH problem 
related to the determinant representation of the correlation functions (jl.3|) and (jl.4|) . 

2.1. Elementary relations 

In this subsection we recall some elementary relations for the correlation functions ()1.3|) 
and (jLil). 

Due to the parity and the time reversal symmetry of the Hamiltonian (jl.lj) the 
functions (jl.3|) and (jl.4|) obey the following symmetry properties 

G h(e) (x,t) = G h{e) (-x,t) (2.1) 

and 

G h{e) (x,t) = G* h{e) (x,-t). (2.2) 

Here the asterisk stands for the complex conjugation. Furthermore, the following 
relation between Gh and G e can be written 

G h (x,t) = -Gl(x,t) + ({^{(x,0^t(0>0)}} (2-3) 
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where {, } stands for the anticommutator. At t — this relation becomes 

G h (x,0) = -G* e (x,0) + 5(x) (2.4) 

which is obvious from equation (jl.2j) . Note also the scaling relation 

G e ( h )(x,t;ii) = y/jlG e ( h )(y/jlx,fit;l) (2.5) 

where G e (h)(x,t; fi) is the correlation function taken at a given chemical potential /i. 

In the following we will assume that x > 0, t > and \i = 1. The correlation 
functions for all other cases will follow from the relations (|2.1j) . ()2.2j) and (|2.5j) . 

U.U. Determinant representation for Gh and G e 

In this subsection, following we write down the determinant representation for the 
correlation functions ()1.3|) and (jl.4|) . We use the notations of j^j. 

In paper |2] the following representation of the correlation functions (jl.3|) and (jl.4j) 
was derived 

G h {x,t) = e —[ d V F(ri)B_4x,t,rj)det(i + V){x,t,r]) (2.6) 

and 

G e (x,t) = -— d V Kl} b ++ (x,t, V )det(I + V)(x,t, V ). (2.7) 

Z7T J_ n 1 — COS 7/ 

To shorten notations we will henceforth omit x, t and r] in the arguments of all functions 
unless they are necessary. 

The objects entering the representations (|2.fi|) and (|2.7jl are defined as follows. The 
function F(rj) is 



c 



= 1 + + -. (2.8) 



The determinant 

V"(/h./ m ) ••• r(A- L ./,- A ; 

dfciv 

V(fcjv,fci) ••• V(fcjv, fcjv) 



jV=0 " •'-J- •'-J- 



;2.9) 



is the Fredholm determinant of a linear integral operator with the kernel 

K — p 

defined on [—1, 1] x [—1, 1]. The functions entering this kernel are defined as follows: 
e+(k) = -^— e - r{k)/2 \{1 - cos r])e TW E Q {k) + sin 77] 

V n (2.11) 
r(k) = ik 2 t — ikx 

roc e -r(p) 

E (/c) = p.v. / dp 



00 



7f(p — fc) 
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Note that V(k,p) is symmetric, V(k,p) = V(p, k), and nonsingular at k — p. 
To define B and b ++ introduce the resolvent operator R : 

I-R= (I + V)-\ (2.12) 

Due to the specific form ()2.10|) of the kernel V(k,p) the resolvent kernel R(k,p) may be 
written as 

R(k>p)= umMpiMm (2 . 13) 

k — p 

where functions f± are the solutions to the integral equations 

f±{k) + J ^ dpV(k,p)f±(p) = e ± {k). (2.14) 
For the pairs e±(k) and f±(k) we will often use the vector notation: 

Now define the functions Bab and C ab 

B ab = J dke a {k)f b {k) C ab = J X dkke a {k)h{k) (2.16) 
where a and b run through two values: a, b = ±. In particular, 

= y dfce_ (*)/_(*;). (2.17) 

The function 6 ++ in (12. 7|) is 

6 ++ = £ ++ - (1 - cosr])G (2.18) 
where (jo(^>*) is the vacuum Green function 

1 roc ( e- 1 "/ 4 ix 2 /4f - / n 

G (s, t) = ^ / dke-^ = J ^ » (2.19) 

2tt J^oo ^ 5(x) £ = 

5. Deformation of the integration contour in equations \2. 6)) and \2. 7| ) 

In this subsection we examine the analytic properties of the functions entering the 
representations (J2.6|) and (J2.7)) . Using these properties we transform the integrals in 
equations ()2.6|) and (|2.7|) to the form convenient for the asymptotical analysis. 
Let us change the integration variable rj in ()2.6j) and ()2.7|) to 

z = e iri . (2.20) 

Then 

e _it C dz ^ 

G h (x,t) = — & —F(z)B__(z)det(I + V)(z) (2.21) 
8vri y w=1 z 

e il f dz - 

G e (x,t) = - i ? -^)& ++ (z)det(/ + \/)(z) (2.22) 

7T1 J\ z \ =l (Z - iy 
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where 



F(z) = 1 + 



+ 



1 



z 2z-l 

1, is oriented counterclockwise. 



(2.23) 



(2.24) 



The integration contour, \z\ 
Introduce the function 

f°° d»e" r(p) 

Eik)= .L£—k ImM0 

and denote by E + (k) and E_(k) its analytic branches in the upper and lower half-planes, 
respectively. The analytic continuations of E + (k) and E_(k) to the real axis satisfy 

E + {k) + £L(Jfe) = £b(*0 £+0) - E-{k) = ie" r(fc) feel. (2.25) 

Using ()2.25|) the kernel V(k,p) = V(k,p; z) can be rewritten as 

V{k,p; z) = {z- l)[Wi(k,p) + z~ 1 W 2 {k,p)] (2.26) 



where 



W 2 (k,p) 



2ty 



exp 



r(k) + r(p) 



2tt 



exp 



r(k) + t(p) 



E + (k)-E + jp) 
k — p 
E_(k)-E_(p) 
k — p 



(2.27) 
(2.28) 



By using the Fredholm theory of linear integral equations (see, e.g., [H], chapter 3, section 
7) the following two facts may be shown: (a) the Fredholm determinant det(J + V)(z) 
with the kernel V(k,p;z) defined by ()2.26)) - (j2.28j) is analytic in the complex z-plane 
except at most at z = and z = oo. (b) the corresponding resolvent kernel R(k,p;z) 
may be represented as 

D(k,p; z) 



R{k,p; z) 



1)- 



(2.29) 



det(I + V)(z) 

where the function D(k,p; z) is analytic in the complex z-plane except at most at z = 
and z = oo. The explicit form of D(k,p; z) is given in the reference cited above. We will 
not use it in our calculations. 

Consider equation (|2.21|) . The function B , entering this equation, is defined by 

equation (j2.17j) . Inverting (|2.14|) 

■I 



e_(*) 



dpR(p,k;z)e-(p) = f-(k) 



and substituting the resulting expression in ([2. 17)1 . one gets 



B__(z) = J dke 2 _(k)-J dk j 



dpe_{k)R{k,p; z)e_(p). 



(2.30) 



(2.31) 



Combining equations ()2.29|) and (|2.31|) one sees that B (z) det(J + V)(z) is analytic in 

the complex ^-plane except at most at z = and z = oo. Keeping this in mind, deform 
the integration contour in (|2.21|) to a contour 7 shown in figure ^ Since the function 
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Figure 1. Deformation of the integration contour. The dashed circle shows the 
integration contour in equations 1)2.211) and (|2.22|l : the solid circle of radius r < 1/2 
shows the integration contour 7 in equations (|2.32l) and (|2.33() . Both contours are 
oriented counterclockwise. The cross indicates the pole of the function F(z) at z = 1/2. 



F(z) has a simple pole at the point z = 1/2, one gets 

G h (x, t) = ^S__(l/2) det(/ + t>)(l/2) 

e _i * f dz - 
+ — i — F(z)B — (z) det(J + V)(z). (2.32) 

07T1 J 7 Z 

Next consider equation (|2.22|) . Like the function B (z) det(/ + V)(z), the function 

b ++ (z) det(I + V)(z) is analytic in the complex z-plane except at most at z = and 
z = 00. Note that despite the singular factor [z — 1)~ 2 the integrand does not have a 
singularity at z — 1 due to the second order zero of b ++ (z) at this point. Indeed, one 
can see from (j2.11|) that e + ~ rj for rj ~ 0. Therefore, ~ 77 for 77 ~ 0, which is seen 
from equation (j2.14|) . Using the definitions ()2.1fij) . ()2.18j) and ()2.2()|1 one can see that 
b ++ ~ (1 — z) 2 for 1 — z ~ 0. Deforming the integration contour in ()2.22j) to the contour 
7 shown in figure ^ one gets 

G e (x,t) = 4e i *6 ++ (l/2) det(J + V)(l/2) 

e il f dz - 
+ -f 7 ^F(z)b ++ (z) det(J + U)(z). (2.33) 

7T1 J 7 (Z - l) 2 

Consider a special case of t = 0. The kernel (j2.28J) becomes equal to zero since 
E_{k) = 0. Hence, the kernel f!2.26|) becomes a linear function of z. This implies that 
(see, e.g., [6 , chapter 3, section 7) (a) the determinant det(J + V)(z) becomes an 
analytic function of z for all z G C (b) the function D(k,p;z) in equation ()2.12|) also 
becomes an analytic function of z for all z G C. Thus, by virtue of 1)2.31)1 the function 

B (z)det(J + V)(z) is analytic in the complex z-plane and the second term in the 

right hand side of ()2.32)) vanishes: 

G h {x,0) = j5__(l/2)det(/+y)(l/2). (2.34) 
The function G e (x, 0) can be obtained from equation ()2.34)1 using equation ()2.4)) . 
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2.4- Differential equations 

In this subsection we write differential equations for the Fredholm determinant det(J + 

V)(x,t,rj) entering the representations (j2.fi j) and ()2.7j) . 

First, we obtain the differential equation with respect to 77. As follows from ()2.20|) 

and ()2.26)) - (j2.28|) the rj derivative of the kernel (j2.1()j) is given by 

i 1 _ e -«? 

d v V(k,p) = -V(k,p) + : W 2 (k,p). (2.35) 

1 — e~ ir ? 1 

Differentiating the identity 

In det(/ + V) = ti ln(7 + V) (2.36) 
we immediately get 

i 1 _ e -"7 „ 
d n lndet(I + V){x,t,7)) = tr[(J + V)~ V] + tr[(/ + V) W 2 ]. (2.37) 

Transform the right hand side of ()2.37|) employing (j2.12j) . (|2.13|l and (|2.15jl . The result 
for the first term is 

tr[(J + V)- 1 ]/] =trR = -i j &k[f(k)} T o- 2 d k f(k) (2.38) 
where a 2 is the second Pauli matrix 

-(:;) -(?o) 

The second term in the right hand side of ()2.37j) is more complicated: 
tr[(J + Vy l W 2 \ = tr[(J - R)W 2 ] 

= J dkW 2 (k,k)-J dk J d P R(k,p)W 2 (p,k). (2.40) 

The differential equations for det(J + V)(x,t,rj) with respect to x and t can be 
obtained in a similar way. They read 

d x lndet{I + V)(x,t,r]) = LB+_ = LB_+ (2.41) 
d t lndet(J + V)(x, t, T]) = -i + + (1 - cosr])G B__] . (2.42) 

Finally, the following initial conditions for the differential equations follow from 
(irmj) and fl2~TTj) : 

det(I + V)(x,t,0) = 1 (2.43) 

and 

det(J + VO(0, 0, 77) = 1. (2.44) 

The differential equations ()2.37|) . (|2.41|) and (|2.42j) will be used in the subsequent 
asymptotical calculation of det(J + V). 
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2.5. Riemann-Hilbert problem 

In this subsection we briefly recall the notion of the Riemann-Hilbert problem (RH 
problem) and its relation to the objects entering the representations (j2.fi j) and (|2.7j) . 

The standard formulation of the Riemann-Hilbert problem is as follows |H [7] : for 
a given oriented contour E in the complex fc-plane and for a given matrix x(k) defined 
on this contour, seek a matrix Y(k) such that 

(a) Y(k) is analytic in C \ E 

(b) Y+(k) = Y_(k) X (k) fceS (2.45) 

(c) Y{k) — > / as A; — ► oo 

The contour E is called 'conjugation contour' or 'jump contour'; the matrix xfo) * s 
called 'conjugation matrix' or 'jump matrix'. The matrices Y±(k) are the boundary 
values of Y{k) as k approaches E : 

Y±(k) = lim Y(k') where k' 6 ± side of E. (2.46) 

k'—>k 

By convention, the + side (respectively, — side) of E lies to the left (respectively, right) 
as one traverses the contour in the direction of the orientation. An example of such an 
oriented contour is shown in figure El 

By using the Cauchy transform one can see that the RH problem 1)2.45)1 is equivalent 
to the matrix- valued integral equation ([7j, chapter XV) 

(The functions Y±(s), x( s ) an d Y±(s)x(s) are supposed to be integrable along the 
contour E. We will impose this condition for all RH problems appearing below). Suppose 
the jump matrix x(k) depends on the external parameter, x, and 

X(k; x) = I + 0{x~ a ) a>0 as x -> oo (2.48) 

uniformly for fceE. Then, it can be shown solving (J2.47)) perturbatively that 

Y = I + 0(x~ a ) as x -> oo (2.49) 

uniformly outside an arbitrarily small vicinity of the contour E. We will use the estimate 
1)2.49)1 in the subsequent asymptotical analysis. 

Next we relate the RH problem ()2.45)) to the objects entering the determinant 
representations ([2.6)1 and ()2.7|) . Consider the following conjugation matrix 




-e-(k)e+(k) e\{k) 
-e 2 _{k) e + {k)e.{k) 



Mk-,x,t)=[ ; y ]-27ri( I (2.50) 



or in the vector notation ()2.15)1 

v {k;x,t) = I + 2ire(k)[e(k)} T o- 2 . (2.51) 

Choose the conjugation contour to be the interval [—1, 1] oriented as shown in figure 121 
and denote this contour as Eq. Let Y be a solution to the RH problem 
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Figure 2. The conjugation contour E - The thick line connecting the points k = — 1 
and k = 1 shows the conjugation contour Eo for the RH problem (|2.52|) . The arrow 
shows the direction of the orientation of this contour. 



(a) Y (k) is analytic in C \ Eo 

(b) Y + (k) = Y„(k)v (k;x,t) fceS . (2.52) 

(c) Y(k) — > / as k — > oo 
Then the components / + and /_ of the vector 

f(k) = Y_(k)e(k) = Y + (k)e(k) fceE (2.53) 

satisfy equations (|2.14j) . Indeed, consider equation (J2.47)) where the conjugation matrix 
X is chosen to be the matrix z/ given in equation (j2.51|) . Multiply equation (|2.47|) by 
the vector e(k) from the right. Equations (j2.14j) follow immediately. 

Combining equations ()2.53j) and (|2.47jl where the matrix x is chosen to be v§ write 

Y(k) =( 1 °) + f JP- ( C "(P)/+(P) ~ e +^fM ) (25A) 
{ ] \ 1 J Lv-k \ e-(p)/_(p) ~e + (p)f-(p) )' 1 ' J 

Comparing equations ()2.16|) and (|2.54j) one sees that the functions B ab and C ab can be 
found from the large k expansion of Y(k) 



Y(k)= t +r "o o ++ +i ,7 + X ++ +f n ■ (2-85) 





Consider an arbitrary compact domain X> in the complex plane of the parameter 
z (|2.2U|) . not containing the point 2 = 0. Since the Fredholm determinant in equations 
(12.21)1 and ()2.22j) is analytic in C\{0}, it only has a finite set of zeroes {z n } = {zi, . . . , z n } 
in D. The integral equations ()2.14|) are solved uniquely for any z ^ {z n } in T>. The 
solutions f± of the integral equations ()2.14)1 can be found from a solution Y(k) of the 
RH problem ()2.52|) by using equation ()2.53|) . The functions B a b and C a b can be found 
from the expansion (|2.55|) . The Fredholm determinant det(J + V)(z) is completely 
defined by the differential equations (|2.41|) and ()2.42|) with the initial condition ()2.44|) 
or, alternatively, by the differential equation ()2.37j) with the initial condition ()2.43jl . 
Since the functions B^(z) det(f + V){z) in (l2~2T| and b ++ {z) det(J + V)(z) in (12~2H) 
are analytic in C \ {0}, their values at the points be found by analytic 

continuation. 
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3. Riemann-Hilbert problem at t = 

In this section we solve the RH problem (|2.52jl asymptotically for x —>■ oo and t = 0. As 
it is obvious from ([2.34)1 . to calculate Gh{x, 0) it is sufficient to solve this RH problem at 
a specific value z = 1/2 of the parameter z. Nevertheless, we will construct the solution 
for all complex z outside an arbitrarily small vicinity of the point z = 0. This general 
result will prove useful in the subsequent analysis of the t ^ case. 

In subsection 13.11 we formulate explicitly the RH problem (|2.52jl for t — 0, see 
equation (j3.3|) . In subsections 13 . 21 and !3 .31 following along the lines of jHE], we transform 
the RH problem ()3.3)1 so as to make it convenient for the large x analysis. In subsection 
13.41 we give an explicit solution to the jump relation (|3.3l b) in a vicinity of the end 
points k = ±1 of the contour E . In subsection 13. 51 we construct the solution to the RH 
problem ()3.3j) in the large x limit by matching the solution given in subsection 13.41 with 
the condition (|3.3l c). In subsection 13.61 we summarize the results of this section. 



3.1. Riemann-Hilbert problem at t = 

In this subsection we give an explicit formulation of the RH problem ()2.52|) at t = 0. 
By taking the limit t — > and assuming that x > in (|2.11|) we get for e±(k) 

e+(k) = l -^=e lkx/2 {\ - e ir? ) e_(Jfe) = -^e~ ikx/2 . (3.1) 

2 \/7T \/7l 



The conjugation matrix ()2.5Uj) becomes 



3 ir? _2ie ifca; e ir ' sin 2 



n 



"« a ' : -''- 0) i 2ie-^ 2-e^ " 1 ' (3 ' 2) 

The oriented contour E shown in figure El remains unchanged. Thus, we have defined 
the RH problem for t = 

(a) Y(k) is analytic in C \ E 

(b) Y+(k) = Y_(k)v (k; Xl 0) fceS (3.3) 

(c) Y(k) — > / as k —>■ oo. 
Introduce a matrix 

e 1 ^ 2 sin 2 
l 

2 

Note that P 2 = I. The jump matrix (|3.2|) and the functions (|3.1|1 satisfy the involutions 

PvoWP- 1 = IM-kT 1 (3.5) 
and (note (j2~T5j) l 

Pe(A;) = e(-A;) (3.6) 
respectively. As a consequence, we have the involution 

PYi-^p- 1 = Y(k) (3.7) 
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Figure 3. The lens-shaped contour Ei ens . The thick lines show the lens-shaped contour 
Siens = ^i U S 2 U £3. It is the conjugation contour for the RH problem This 
contour divides the complex fc-plane into three regions: Oy, and Sloo- 



for the solution Y{k) of the RH problem ()3.3|) and, by virtue of ()2.53|) . the involution 
Pf(k) = f(-k). (3.8) 

for /(*). 

3.2. Riemann- Hilbert problem on a lens-shaped contour 

In this subsection an RH problem equivalent to the RH problem ()3.3|) is formulated on 
the lens-shaped conjugation contour shown in figure El 
Represent Y(k) product of two matrices 

Y(k) = U{k)L{k) (3.9) 

and define L(k) by the formula 

{I k G fioo 

u x (k) ken v . (3.10) 

A lens-shaped contour Si ens = S1UE2UE3 dividing the complex fc-plane into the regions 
flu, and Qoo is displayed in figure El The matrices v\ and u 3 are defined as 

1 -2ie ite sin 2 2 \ 

vi(k;x) = I Q 1 I (3.11) 

■*<*;*)=! 2ie -i e ->, J)- (312) 

Introducing a matrix u 2 by the formula 

■*= (fe-°») p- 13 ' 

one can see that the following representation for the conjugation matrix ()3.2|) can be 
written 

^0 = ^2^1- (3.14) 



Zero temperature correlation functions for the impenetrable fermion gas 



12 



Note that 

P Vl (k)p- 1 = [^(-k)]- 1 Pv 2 P- 1 = v 2 l (3.15) 

where the matrix P is defined by ()3.4j) . 

It follows from equations ([3~3]l . (Q , (EHTH) and (I3~T1) that £/ (fc) solves the following 
RH problem 

(a) U (k) is analytic in C \ Si cns 

(b) U+(k) = U-(k)ui(k) ke^i i = 1,2,3 (3.16) 

(c) U{k) — > / as k —* oo 

We will calculate the large x asymptotics of £/(&); the asymptotics of will then 
follow from (|3.9|) . The advantage of the RH problem ()3.16|) from the point of view of 
the large x analysis will be clear from the next subsection. 

3.3. Factorization of the Riemann-Hilbert problem A3. 16}) 

In this subsection we present the scheme of the analysis of the RH problem ()3.16|) in 
the large x limit. 

The x-dependent entry of the jump matrix V\ (respectively, z/ 3 ) entering (|3.16|) 
decays exponentially as x — > oo everywhere on the jump contour Si (respectively, S3) 
except at the points k = ±1. The jump matrix z/ 2 is a constant matrix. For the RH 
problems with such a behaviour of the jump matrices, the following scheme of the large 
x analysis can be employed [HE]- Denote by V a vicinity of the points k = ±1. Split 
U(k) in a product of three matrices 

U = KSQ. (3.17) 

The matrix Q(k) satisfies the jump relation Q+(k) = Q-{k)i>2 everywhere on the contour 
S 2 and the jump relation ([3. 161 b) on the part of contour Si ens lying in V. The matrix S(k) 
represents the corrections to Q(k) due to the mismatch between the analytic branches 
of Q(k) in and outside V. The matrix K(k) represents the corrections to S(k)Q(k) due 
to the discontinuities, exponentially small as x — > 00, across the parts of the contours 
Si and S 3 lying outside V. Below we give the precise sense to Q, S and K. 

Let V consist of two discs T>r and T>l shown in figure 0] Let the radii of these 
disks be equal and less than one. Denote the domain C \ (T>r UVl) as V^,. Introduce 
a matrix Q(k) by 

( Q R (k) keV R 

Q{k) = I Q L {k) kEV L (3.18) 

[ Qoz(k) kev^ 

where the matrices Qoo, Qr and Ql are defined as follows. The matrix Qoo(k) is analytic 
everywhere in T)^ except the contour S 2 H V^. On this contour it satisfies the jump 
relation with the jump matrix (|3.13J) 

[Qoo(k)] + = [Qoo(A;)]_z/ 2 k E S 2 n X?oo (3.19) 



Zero temperature correlation functions for the impenetrable fermion gas 



13 




and is normalized at infinity: 

Qoo(k) -> I as k -> oa. (3.20) 
We choose the following (obviously, not unique) solution to (|3.19|) and (j3.2U)l : 

Note that it satisfies the involution 

Q 00 (k) = PQ 00 (-k)P- 1 . (3.22) 

The matrix Qn(k) is analytic everywhere in Dr except the contour Si ens n Dr. On this 
contour it satisfies the jump relation with the jump matrices ()3.11|) - ()3.13|) 

[Q R (k)] + = [Q R (k)]_u i (k) ke^nVR z = l,2,3. (3.23) 

A solution to the jump relation ()3.23|) is given in subsection 13.41 The matrix is 
analytic everywhere in D £ , except the contour Si ons HD L . On this contour it satisfies 
the jump relation similar to (J3.23|) 

[QL(k)] + = [Q L (k)]_Ui(k) keZi(lV L z = l,2,3. (3.24) 

A solution to the jump relation ()3.24j) is given in subsection 13.41 

It should be stressed that any solution to the jump relation (|3.23|) (respectively, 
(13.24)0 multiplied from the left by an arbitrary matrix with the entries being entire 
functions of k in D R (respectively, V L ), still remains a solution to (|3.23|) (respectively, 
(13.24)1 ). In subsequent calculations we will employ this property to simplify the form of 
the conjugation matrices Or and 6 L defined by equation (J3.27)) below. 

Next define the matrix S entering the decomposition ()3.17j) as a solution to the 
following RH problem: 

(a) S(k) is analytic in C \ (dV R U dD L ) 

(b) S + (k) = S-(k)0R >L (k) kedDR L (3.25) 

(c) S(k) — > / as k — >• oo 

with the clockwise oriented conjugation contours &Dr^ shown in figure |U The jump 
matrices 9R^(k) are chosen to ensure the continuity of SQ across &Dr^l '■ 

S+(k)Q + (k) = S_(k)Q_(k) kedV R>L . (3.26) 
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£1 




4> 1 




S 3 



Figure 5. The jump contour for the jump relation l|3.23[) . The disk T>r is shaded 
in gray. The segments of the contours Si and S3 lying in T>r are straight and 
perpendicular to S 2 . 



This implies 



R ,L(k) = Q R , L (k)[Q oo {k)] 



k e dV 



•r,l- 



(3.27) 



We solve the RH problem ()3.25|) in the large x limit in subsection 13.51 

With the matrices Q and S defined above one can easily see that the matrix K in 
the decomposition (J3.17|) solves the following RH problem 

(a) K(k) is analytic in C \ (£1,3 fl £>oo) 

(b) K+{k) = K^{k)[S(k)Q(k)]u l (k)[S(k)Q{k)}- 1 fceE^Doo 2 = 1,3 (3.28) 

(c) K(k) — > / as k —>■ 00 

We solve the RH problem ()3.28|) in the large x limit in subsection 13.61 
3.4- Exact solutions to the jump relations \3. 23}) and \3. 24 ) 

In this subsection we explicitly construct matrices Qr and Ql satisfying the jump 
relations (|3.23J) and ()3.24|) . respectively. 

Since the conjugation matrix ()2.50|) is 27r-periodic in Kerj we will restrict our 
considerations to the interval 



Consider the part of the contour Si ens lying in the disk T)r. For the convenience 
of presentation we choose Si and S 3 to be perpendicular to S 2 , as shown in figure El 
Consider a matrix $#(£;; x) with the entries 



7r < Ker/ < 7r. 



(3.29) 



To shorten notations, introduce the parameter 
_ t |Rer/| 



(3.30) 



7T 
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®u(k;x) 



x~^si(k - l)eT* — 1; -ix(k - 1) 



_n_ nqUk — l)e 4 e 

-T 2tv : - 

P(£)sin 2 f 

3i2 



-ix(k — 1) 



__n_irs 2 (k - l)e * e lx 

X 2 - _„ ; „ ; * 



(3.31) 



r 2 (-&) 



l + -i-,l;ix(A;-l) 

Z7T 



x^(, 2 {k - lje^f* — , 1; iar(A; — 1) 

-27T 



Here ^(a, b;w) is the Tricomi \l/ function, which has a branch cut discontinuity across 
the negative real axis [B] 

(a, 1, tw + iO) - [1 - 2ie~ ia7r sin(Tra)] *(a, 1, w - iO) 
27rie" iaw 



r 2 fa) 



-e M *(l - a, 1, -tu) we 



and Si,2( w ) are piecewise-constant functions defined by 



and 



-7r < argw < — g 
-~ < argu> < 7r 



— 7r < argu> < | 
~ < argw < 7r 



(3.32) 



(3.33) 



(3.34) 



Using ()3.32|) it is easy to check that $ R (k;x) satisfies the jump relation (j3.23j) on the 
contour Ei cns R V R shown in figure El 

[$ R (k',x)] + = [$ R (k;x)]-i/i(k;x) ke^UV R i = 1,2,3. (3.35) 

It follows from the involution ()3.15|) that a solution $l(&) to the jump relation ()3.24|) 
may be obtained from <& R (k) via the involution 

$ L (k; x) = P$ R (-k; x)P- 1 k E V L . (3.36) 

As was mentioned in the paragraph following equation (|3.24j) . any matrix Q R of 
the form 

Q R (k) = E(k)$ R (k) (3.37) 

with E{k) being an arbitrary analytic in T> R matrix, satisfies ()3.23j) . Choosing 

(Jfc + l)-& 
(Jfe+1)& 

one ensures that Q R {k) goes to Qoo(k), as x — > oo, for any e ^ 0, see f!3.41jl . This will 
facilitate the asymptotic solution of the RH problem ()3.25|) . Taking into account (|3.36|) 
and ()3.37|) define Qhik) in T>l by the formula 

Q L (k) = PQ R (-k)p- 1 k e V L . (3.39) 



E(k) 



(3.38) 
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To calculate the large x asymptotics of Q R (k) we use the asymptotic formula for 
the Tricomi \1/ function (|H|, section 6.13.1) 



Applying fl35PD to (EPT1) one gets from (EP7I) : 





f" a 2 




w~ a 


1 






w 





as io — > oo. (3.40) 



where 



-i2 - 

i7re 2e 



2 



sm 



The large x asymptotics of Qhik) follows from (J3.39)) and (|3.41jl . Note that the terms 
in asymptotic expansions of Q R explicitly written down in (|3.41j) do not have jumps 
across the contours Ex and E 3 . This reflects the fact that the x-dependent entries in 
the jump matrices v\ and u 3 are exponentially small on Ei and E 3 , respectively. 

3.5. Solution to the Riemann-Hilbert problem \3. 25\) in the large x limit 

In this subsection we solve the RH problem ()3.25j) in the large x limit. 

Consider the conjugation matrices 9 R and 9l of the RH problem 1)3.25)1 . They are 
defined by the formula ()3.27|) and satisfy the involution 

L (k) = P9 R (-k)p- 1 k G dV L (3.43) 

as can be shown using the involutions (J3.22)) and (|3.39|) . Substituting the asymptotic 
expression (j3.41|) into ()3.27|) one gets the approximation of the matrix 9 R by a matrix 
Or 

6 R =~e R [I + 0{x- 1 -*)) (3.44) 



uniform for k G &D R . Here 

x{k - 1) 

and 



9 R (k) =1+ V f f k \, k G &D R (3.45) 



b R [x(k + l)]« J 
It follows from the involution ()3.43|) that the matrix 9l is approximated by a matrix 9l 
9 L = 9 L [I + 0(x- 1 - e )} (3.47) 
uniformly for k G &Dl, where 

9 L {k) = P9 R (-k)P~ l k G dV L . (3.48) 
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Consider the RH problem 

(a) S(k) is analytic in C \ (dV R U dV L ) 

(b) S + {k) = S4k)9 R>L {k) ke&D R>L (3.49) 

(c) S(k) — > I as k — ► oo 

with the matrices # R and 6 L defined by equations ()3.45|) and ()3.48|h respectively. 
Comparing (jHHj) and (pPTjl with (J213J) and (J21T?|) we conclude that the solution S(k) 
of the RH problem ()3.49|) approximates the solution S(k) of the RH problem ()3.25|) 

S = S[I + Oix- 1 - 6 )] (3.50) 

uniformly outside a vicinity of the conjugation contours dT> R and 5X>£. 

Due to a simple analytical structure of the conjugation matrices (I3.45j) and ()3.48j) 
the RH problem ()3.49|) can be solved exactly. To begin with, introduce the following 
notation for the analytic branches of S(k) 

( S R (k) keV R 
§(k) = I S L {k) keV L . (3.51) 
[ Socik) keV^ 

The matrix 8 R (k), analytically continued from dT> R into T> R has a simple pole at k — 1, 
as follows from the representation ()3.45|) . The matrix S R (k) is analytic in T> R . Therefore, 
the jump relation ()3.49l b) 

~Soo{k) = S R (k)6 R (k) k G dV R (3.52) 

implies that the analytic continuation of Soo(k) into T> R has a simple pole at k — 1. 
Similarly, the analytic continuation of S^k) into T>l has a simple pole at k — — 1. 
Finally, taking into account the condition given in equation (I3.49l c). S^oo) = I, we 
see that the most general form of Soo(k) is 

~S 00 {k) = I + + keC (3.53) 
k — 1 k + 1 

with the matrices A R and being independent of k. 

The involution 

PS^P- 1 = S^-k) (3.54) 
which follows from the involution (J3.48|) . provides a simple relation between A R and Ai 

A R = -PAlP- 1 (3.55) 
so we only need to calculate one of the matrices, say, A R . For this, invert (|3.52j) : 

~S R (k)=~S OQ {k)i6 R {k)Y 1 kedV R . (3.56) 

Using ()3.45|) and (|3.53j) the right hand side of equation (|3.56J) can be analytically 
continued into V R . In the following, equation (|3.56|) will be understood to be such 
a continuation. The matrix [^(fc)] -1 has the form 

vn(k) 
x(k — 1 



[^(^-^[detMfc)]- 1 



(3.57) 
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where 



dBtgH(AO = l + ^* )a (3.58) 
and det vr is a fc-independent number 

det "«=(s) 4 + (s) 2 - (3 - 59) 

The right hand side of ()3.56|) must be analytic in since the left hand side is. The 
matrix [^(A;)]" 1 goes to zero as k — > 1, cancelling the pole of S^^k) at this point. 
However, is singular at the points where 

det 9 R (k) = 0. (3.60) 

Equation ()3.6()|) is quadratic in k. Its solutions, /cf 2 , lie inside the disk Vr for sufficiently 
large x. Therefore, to cancel the singularities in ()3.56|) the matrix Ar has to satisfy the 
pair of matrix equations 



A R P^ArP 

' UR _ 1 t,R I i 
«,j 2 1 "-1,2 ^ - 1 



0. (3.61) 



The determinant of the matrix in square brackets in (|3.61|) is equal to det6^(fcf 2 ) 
and, therefore, is equal to zero. Employing this fact one gets (after some algebra) the 
matrix Ar. The large x asymptotics of Ar is 

Ar = l( * - Uif -W» si, l ^ + , 



where 



and 



Tre^x)" 1 "- 7re- i:E (2a;)~ 1+ - 

r (-^) S]n l r (If) 8111 I 



K = 1 + k\ + k\. (3.64) 
Applying the involution ()3.55j) to ()3.62|1 one gets the large x asymptotics of Al : 

: 4 :,F ( *f '?f r;'L W-'>- (3.65) 




5. 6. Riemann- Hilbert problem |X3J).' results in the large x limit 

In this subsection we show that the the matrix K in the decomposition ()3.17|) does not 
contribute to the asymptotic solution of the RH problem (|3.16|) to the order we are 
interested in. This will conclude the asymptotic solution of the RH problem ()3.3|) . 

Consider the decomposition (|3.17j) . In subsection 13.51 we constructed explicitly 
the approximation S to S up to the order of x~ 1 ~ e . This approximation is uniform 
outside an arbitrarily small vicinity of the conjugation contours &Dr and &D^. The 
matrix Q was constructed explicitly in subsections 13.31 and 13.41 Using these results, 
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solve asymptotically the RH problem (|3.28|) for the matrix K. The matrices V\ and z/ 3 , 
entering ()3.28l b). converge as x — ► oo to the identity matrix uniformly and exponentially 
fast on the contours Si n T)^ and E3 D £>oo, respectively. One can prove that 

K = I + o(x'°°) (3.66) 

uniformly outside a vicinity of the contours S 13 D Pqo- The basic idea in obtaining 
()3.66|1 is the same as in obtaining (|3.5()j) : to use (|2.47jl . One should take into account, 
though, that we have now, contrary to (|3.50j) . no uniform estimate for the jump matrix 
in ()3.281 b) in a vicinity of the contours &Dr,l since we have no uniform estimate for 
S(k) entering this matrix, in this vicinity. For further details we refer the reader to 
reference jl]. 

Substituting ()3.5()|) and ()3.66|) into the decomposition ()3.17j) one gets the estimate 

U = SQ[I + 0(x- l - e )} (3.67) 

uniform outside an arbitrarily small vicinity of the contours dV^L and £1,3 fl Poo- Using 
(13. 9|) one gets the approximation to Y{k) up to the order of x~ 1_e , uniform outside this 
vicinity. This completes the solution to the RH problem ()3.3|) in the large x limit. 



4. Equal-time correlation functions 

In this section we calculate the equal-time correlation functions Gh{x, 0) and G e (x, 0) 
in the large x limit. The asymptotic expressions for these functions were obtained in 
papers [211101 ■ The main purpose of this section is a to give a detailed illustration of 
our method on this known example. The procedures of this section will be extensively 
used in the subsequent analysis of the time-dependent correlation functions. 

In subsection 14.11 we calculate the functions B a b ()2.16)1 from the large k expansion 
()2.55|) of Y(k). From thus obtained B + _ we calculate the Fredholm determinant 
det(J + V) in the large x limit, see equation ()4.6j) . The constant C(rj) entering ()4.6|) is 
calculated in subsection 14.21 We give the answer for Gh{x, 0) and G e (x, 0) in subsection 
IP1 

4-1. Large x asymptotics of B a b and of the Fredholm determinant det(/ + V) 

In this subsection we use the results of section |3] on the asymptotic solution of the RH 
problem ()3.3j) to calculate the functions B a t and the Fredholm determinant det(J + V) 
in the large x limit. 

It follows from (EOTTjl . flUHj) , lETTEl) and (l33T| that in the vicinity of k = 00 the 
following uniform approximation is valid: 

Y = S OD Q OD + O(x- 1 ^). (4.1) 

Recall that the matrix Qoo(k) is given by the formula (|3.21|) . the matrix S^k) by the 
formulas ()3.53j) . ()3.62|) and ()3.65j) . The approximation (j4.1j) being uniform, the functions 
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B ab can be calculated from the large k expansion (|2.55|) of equation (|4.1|) . Expanding 
(13.21}) and ()3.53j) in inverse powers of k and comparing with ()2.55|) one gets 

-_ B B + f + + )-A R + A L - V -a 3 + 0( X -^) (4.2) 

where 03 is the third Pauli matrix ()2.39j) . Substituting equations ()3.62j) and ()3.(i5j) into 
()4.2|) one gets for B + _ and B 



B+_ = i + - 

7T K 



Xn Xr H 

* L X V27T 



\ 2 



+ O (aT 1 " 6 ) (4.3) 



= 2ie~^ esc |(xt - xr)[1 + O (oT e )] (4.4) 

where and are defined by ()3.63|1 and k is defined by (j3.64jl . 

To calculate the Fredholm determinant det(J + V") we integrate (|2.41|) taking into 
account the boundary condition (|2.44|) : 

lndet(J + y) = i f dyB + -(y). (4.5) 
Jo 

Substituting equation ()4.3|) into ()4.5|) and calculating the integral in the large x limit, 
one obtains 

lndet(J + y) = l ^x-2 (JL) 2 ]nx + C(ri) + 0{x~ e ). (4.6) 

7T V27T/ 

We calculate the integration constant C{rj) in subsection 14.21 
4-2. Calculation ofC{rj) 

In this subsection we calculate the constant C(rj) in (|4.fi|) using the differential equation 

Consider equation ()2.37j) . For t = the kernel W2(k,p) vanishes and, therefore, the 
second term in the right hand side of ()2.37|) vanishes. The remaining term is given by 
equation (J2.38J) . With the help of the involution ()3.8|) and the identity 

P T a 2 P = -a 2 (4.7) 

the integral in ()2.38|) can be reduced to the integral over positive k, yielding 

<Vndet(/ + t>)(x,0,77) = I dk[f(k)] T a 2 d k f(k). (4.8) 

1 ~~ e Jo 

The large x asymptotics of the function f(k) in equation ()4.8|) is calculated as 
follows. Use equation ()2.53|) 

f{k) = Y+(k)e(k) fcGE 2 (4.9) 

and replace Y(k) by its large x asymptotics obtained in sectional It follows from (|3.9|) 
that 

Y+(k) = U + {k)v x {k) fceE 2 (4.10) 
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Figure 6. Analytic continuation. Equation (|3.67|) is analytically continued from the 
disk T)r into the region Q,r containing the interval [0, 1]. 



where U(k) solves the RH problem (j3.16j) . The formula (j3.67j) gives the approximation to 
U{k) uniform outside an arbitrarily small vicinity of the contours dV RL and E li3 HV^. 
Since the radius of V R was chosen initially to be less than one, the contour dV R crosses 
the interval [0, 1]. The integration variable k in (I4.8J1 lies in this interval and, therefore, 
the estimate (jH.fi 7j) can not be applied everywhere on the integration contour. This 
problem is circumvented by the analytic continuation of equation (j3.67j) from the disk 
T> R into a bigger region containing the interval [0,1], as shown in figure The 
matrices Qn(k) and S R (k) are given in Q R by the same formulas as in T> R : the matrix 
Qr by (pUTTJl . (l3~3%J) and (pPTl) : the matrix S R by (|3*3fl . Using flOJ) and (ITTUl) one 
gets the approximation to f(k) 

f=~S R g[I + 0{x- 1 -*)} (4.11) 

uniform for k G [0, 1]. The function g is defined by the equation 

9(k) = ( ^ J = [Q R (k)} + Mk)e(k) k G [0, 1]. (4.12) 

The matrix ui entering (j4.12)l is given by equation (jH.llJ) . the function e by (j2.15j) and 
(l3~Tl) . the matrix Q R (k) by (l3~3TJ) . IET37I) and 1CT3HI) . With the help of the identity jS] 

^(a, b; w) = J u {b) ^ aX ^(a, b; w) + ^&^(^ e w ^(b - a, b; -w) (4.13) 
1 (o — a) 1 (a) 

where A = sgn Imu>, one gets for (j4.12j) 



y^e 4 _ L 



g+(k) = - -^— e-^ +ix x-^(k + ly^xFx + 1, 1; Lr(fc - 1 



r (-£) 

y/Tie~ ~i ikx 



V 



2tt 



(4.14) 



g-(fc) = - , = v . v e 2 i 2 *(Hl) 2 -iFi — ,l;ix(k-l 
r(^)smf L2vr 

Using the uniform estimate (j4.11jl and the formulas (j4.14j) one can derive the 
following estimate for (j4.8j) 

d v \ndet(I + V) = / dkfiWfatdkftkftl + Oix- 1 -')]. (4.15) 

1 — e Jo 
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Performing the integral in (J4.15|) asymptotically one gets 



where 



0„]ndet(/ + y) 

CM 



IX 7] 
7T 7T 



V 

2tt 2 



- ±hix + CM + O (x~ e ) 

- 2(ln2 + 1) 



and ip(w) is the digamma function, 

ib(w) = - — lnlYiu). 
aw 
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(4.16) 
(4.17) 

(4.18) 



Integrating equation ()4.16|) over rj, taking into account the initial condition ()2.43|) and 
comparing the resulting expression with (J4.6)) one gets 



C{ri) = -Mn 

7T 



r(i-&) 



ri 2 1 
' (ln2 + l)-- 



2tt 



7T 



dA In 



r(i-A) 



(4.19) 



4-3. Equal-time correlation functions: the results 



In this subsection we give the answer for the correlation functions Gh(x, 0) and G e (x, 0) 
in the large x limit. 

The correlation function Gh(x, 0) is given by the formula (|2.34j) . The asymptotics 

of B is given by (|4.4|) : the asymptotics of det(J + V) by ()4.6|) with the constant C(rj) 

given by ([4. 19)1 . Recall that z = exp(i^), therefore the point z = 1/2 corresponds to 
rj = iln2. The final answer for Gh{x, 0) is 



Gh(x, 0) = Ex "sin ^x — x In x ) exp 

where the anomalous exponent a is given by 

2 



In 2 



-x 



71 



a 



1 



1 /ln2 



2 \ 7T 

the phase shift xo by 

(ln2) 2 



x 



- 21m 



7T 



lnT 



iln2 
"2^T 



and the constant H by 

S = -47^ exp <j C(iln2) - 2Re 



lnT 



iln2 
~2t7 



(4.20) 



(4.21) 



(4.22) 



(4.23) 



where C(r]) is defined by equation ([4.19)1 . The relative correction to (|4.20|) is of the 
order of x~ x . Note that equations ()4.20j) and (|4.21j) were first derived in paper jHj. The 
phase shift (14.22)1 and the constant ()4.23)1 were calculated in paper [TU] . 

The asymptotics of the function G e (x, 0) is obtained from equation ()4.20jl by using 
equation (|2.4|) . Recall that in the derivation of equation ()4.20jl we assumed that x is 
positive. The result for negative x follows from equation (|2.1)1 . 
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Figure 7. The conjugation contour Eo U S s for the RH problem (|5.10Jl . The sectors 
iSi and Sm are shaded in gray. 

5. Riemann-Hilbert problem at t ^ 0. Space-like region 

In this section we find the asymptotic solution of the Riemann-Hilbert problem ()2.52|) 
for |x| > 2\t\. We assume that x > and t > 0. We will consider the large x asymptotics 
for a fixed value of 

In subsection 15.11 we reformulate the RH problem (|2.52j) so as to make it similar 
to the RH problem ()3.3|) . In subsections 15.21 and 15.31 we repeat the construction made 
for the t = case in subsections 13.21 and 13. 3^ respectively. In subsection 15.41 we give an 
explicit solution to the jump relation ([5 .101 b) in a vicinity of the end points k = ±1 of 
the contour E . In subsection 15 . 51 we perform the matching procedure analogous to that 
of subsection 13.51 In subsection 15.61 we summarize the results of this section. 



5.1. Reformulation of the Riemann-Hilbert problem h2.5ty 

In this subsection we reformulate the RH problem ()2.52|) so as to make it similar to the 
RH problem Q . 

Represent Y{k) clS cL product of two matrices 

Y{k) = Z{k)V{k) (5.2) 

and define V(k) by 

m = fi -2Wi[ B (*)-ie-^ W ]\ (53) 

with E(k) given by equation ()2.24j) and g(k) being a piecewise-constant function 

( 1 k G Si 

g(k) = < -1 k e Sm ■ (5-4) 
j otherwise 

The sectors Si and Sm are shown in figure [7| The oriented contour S s separating the 
sector Si from Su and the sector Sm from S\y is chosen so that (a) it passes through 
the point k s (b) the real part of r{k) is positive for all k G S s except k = k s 

Rer(k) > keT, s \ {k s } (5.5) 
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and 

Rer(fc) — > +00 as k — > 00 G E s . (5-6) 

It is easy to see that V(k) is analytic everywhere in the complex fc-plane except the 
contour E s , where it satisfies the jump relation 

V + (k) = \fiL.(k)]- 1 V4k) fceE s (5.7) 

with 

m ( 1 2ie-^sin 2 f \ 

M*0 = I x J • (5-8) 

Note that the condition (|5.6j) implies 

V(Jfc) -> / as k -> 00. (5.9) 
With V(fe) and £ s defined above, one gets the following RH problem for Z(fc) 

(a) is analytic in C \ (S U E s ) 

(b) Z + (k) = Z-{k)m{k) ke^i i = 0,s (5.10) 

(c) Z(k) — > / as k — > 00 

where the contours £0 an d E s are shown in figure [7| Since is analytic across E s , 
the matrix /x s should satisfy (|5.7|) and, therefore, is given by (|5.8j) . The matrix /x can 
be obtained comparing the jump relations (|2.521 b) and (|5.10l b) 

^ (k) =V(k)u (k;x J t)[V(k)}- 1 /eGEo. (5.11) 

Representing (|2.5Uj) in the form 

*-^(r sin2 f- (ft 0(^e°O(r 2si r +w i (-) 

with E ± given by ()2.24|) and ()2.25j) one gets 

/ e 1 " -2ie- T (*)e^ sin 2 § ' 
= 2ie^) 2-e- (5 ' 13) 



We see that the matrix /xo is different from the matrix uq, defined by (|3.2jl . only in that 
i/cx in the latter is replaced by —r(k) in the former. The matrix /xo obeys the involution 

P/x (&; a:, i)^ -1 = M~*; x, -^T 1 (5.14) 

with the matrix P defined by equation (|3.4|) . This involution is a natural generalization 
of the involution ()3.5jl . 




Figure 8. The conjugation contour for the RH problem l|5.22|l . The contour consists 
of the lens-shaped contour Si en s = £1 U £2 U S3 and the contour E s . 



5.2. Riemann- Hilbert problem on the contour Ei ens U S s 

In this subsection an RH problem equivalent to the RH problem (I5.10J1 is formulated on 
the conjugation contour Si cns U E s . The procedure will resemble that of subsection 13.21 
Represent Z(k) clS cL product of two matrices 

Z(k) = U(k)L(k) (5.15) 

and define L(k) by the formula 

/Xi(Jfe) k e ttu . (5.16) 

[^3(^)1 1 kett D 

The lens-shaped contour S lens = Si U £ 2 U S 3 is the same as in subsection 13.21 It is 
shown in figure |H1 The matrices \i\ and ^ are defined as 

/,(*)=(; -^'-'i) (5,r, 

/«.(*) = ( 2i J % -« I ) ■ ("8) 

For the conjugation matrix ()5.13j) the following representation can be written 

£*o = £*3^2£*i (5.19) 
where the matrix 

» = ( A ) (5 - 20) 

coincides with the matrix 2^ defined by equation ()3.13|) . Note that 

P^k-x^p- 1 = ^(-fcjx.-t)]" 1 Pfjt 2 P~ 1 = fh 1 (5.21) 

where the matrix P is defined by ()3.4|) . 

It follows from equations (l5~TUl) . dSHSj), (l5~T7)J) and (0511 that J7(fe) solves the 
following RH problem 

(a) U{k) is analytic in C \ (Si cns U S s ) 

(b) U+{k) = U-{k)yLi{k) keT^i z = l,2,3,s (5.22) 

(c) U(k) ^ I as k ^ 00 
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Figure 9. The contour used in the factorization of the RH problem l|5.22|l . The disks 
Dri T^l aud are shaded in gray. 



where the contours Sj are shown in figure |H1 We will calculate the large x asymptotics 
of U(k); the asymptotics of Z{k) will then follow from (j5.15j) . 

5.3. Factorization of the Riemann-Hilbert problem 

In this subsection we present the scheme of the analysis of the RH problem (|5.22j) in 
the large x limit. The scheme is similar to the one given in subsection 13.31 

Like the t — case, the jump matrix n\ (respectively, fi 3 ) goes to the identity 
matrix, as x — > oo, everywhere on the contour Si (respectively, E 3 ) except at the points 
k = ±1. The jump matrix /i s goes to the identity matrix everywhere on S s except at 
the point k = k s . Therefore, to construct a uniform approximation to U(k) we use the 
same scheme as in subsection 13.31 Introduce three disks T> R , T>l and T> s as shown in 
figure M Denote the domain C \ (V R U V L U V s ) as V^. Represent the solution U{k) of 
the RH problem ()5.22j) as a product of three matrices 

U = KSQ. (5.23) 

The matrix Q(k) is given by the formula 

f Qn,L{k) k e V RjL 

Q{k) = I Q s (k) kEV s . (5.24) 

[ Qoo(k) k G 

The matrix is given by equation (j3.21|) . The matrix Q s is defined in T> s and satisfies 

[Qs(k)] + = [Q s (k)]_fx s (k) kez s nv s . (5.25) 

We choose the following solution of the jump relation ()5.25|) : 

( 1 2 sin 2 2 f ffe e- T(p) \ 
Q s {k) = Qoo(k) I Q zsm 2h s 2, P - k (5 26) 

where the integral runs along the contour S s shown in figure El The matrices Qr and 
Ql are defined in the disks V R and V L , respectively, and satisfy there the jump relations 

[Q R (k)] + = [Q R (k)}_ f i l (k) keZiHVn i = l,2,3 (5.27) 

and 



[Q L (k)] + = [Q L (k)]_^(k) keZ, t nv L i = l,2,3. (5.28) 



Zero temperature correlation functions for the impenetrable fermion gas 



27 




Figure 10. Conformal map Xr. The disk T>r and its image V R are shaded in gray. 
The conjugation contours Si and S3 are chosen so that their images S^ and S 3 are 
parallel to the imaginary axis. 



We give solutions to the jump relations ()5.27|) and ()5.28|) in subsection 15.41 
The matrix S is defined as a solution of the RH problem 

(a) S(k) is analytic in C \ (dV R U dV L U &D S ) 

(b) S+(k) = S-(k)6 a (k) kedV a a = R,L,s (5.29) 

(c) S(k) — > I as k — ► 00 

with the clockwise oriented conjugation contours &D R l s shown in figure E] and the 
conjugation matrices 6r,l,s defined by 

e a {k) = Q a {k)[Q 00 {k)Y 1 kedV a a = R,L,s. (5.30) 

We solve the RH problem ()5.29|) in the large x limit in subsection 15.51 

The matrix K in the decomposition (J5.23)) is constructed in complete analogy with 
the matrix K considered in subsection 13.31 

5.4- Exact solutions to the jump relations \5.21\) and 

In this subsection we solve explicitly the jump relations ()5.27j) and 1)5.28)1 by mapping 
them on the jump relations (j3.23j) and (|3.24jh solved in subsection 13.41 
Introduce the 'light cone' variables xr and xl 

XR = x — 2t XL = x + 2t. (5.31) 

Consider the jump relation ()5.27j) and make the following conformal map 

X R (k) = 1 - r(fc) .~ T(1) k e V R (5.32) 

of the disk T>r in the complex fc-plane onto a vicinity T>' R of the point = 1 in the 
A/j-plane, as shown in figure ITTH Introducing a /c-independent matrix T by the formula 

/ e~ it/2 \ 

T ={ *») (5 - 33) 

one can easily check that 

m{k)=T- 1 u i [X R {k);x R ]T keV R z = l,2,3. (5.34) 
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Here the matrices /xi, /x 3 and /j 2 are defined by the formulas (|5.17|) . (J5.18)) and (|5.2(J|) . 
respectively; the matrices u±, u 3 and u 2 are defined by the formulas 1)3.11)1 . ()3.12jl and 
(13.13)1 . respectively. Choose the shape of the contours Si and S 3 in the fc-plane so that 
their images S' x and E 3 in the A^-plane are parallel to the imaginary axis, as shown in 
figure HOI Write equation 1)3.35)1 with k replaced by Xr and x replaced by xr : 

[<&r{Xr]Xr)} + = [$R( y \R;x R )]_v i ( y \R;x R ) A fi eS; * = 1,2,3 (5.35) 

where $ R is defined by (13.31)1 . By virtue of ()5.34)1 the jump relation ()5.35)) is equivalent 
to 

{$R[\R(k); x R ]T} + = {$ R [\ R (k); x R ]T}_ m(k) fcGS, i = 1, 2, 3. (5.36) 

One sees that the matrix <& R [\ R (k); xr]T solves the jump relation ()5.27j) in the disk D R . 
Now we are in a position to write down an ansatz for Q R (k) : 

Q R {k) = E R (k)$ R [\ R (k); x R ]T (5.37) 

where 

s a <*) = g.w(MVr* [A R (*)°-i]ft ) T_1 ' < 5 ' 38 > 

With this choice of E R {k) the matrix Q R {k) goes to Qoo(k), as i -> oo, for any e 7^ 0. 

Taking into account the involution ()5.21jl choose the following solution to the jump 
relation (jO^I 

Q L {k-x,t) = PQ R {-k-x,-t)p- 1 keV L (5.39) 
where Qr is given by formulas 1)5.37)1 and 1)5.38)1 . 

5.5. Solution to the Riemann-Hilbert problem \5. 2fJ\) in the large x limit 

In this subsection we solve the RH problem ()5.29)1 in the large x limit. The procedure 
will resemble that of subsection 13.51 

Consider the conjugation matrices 6 R and 6 L defined by equation ()5.30j) . They 
satisfy the involution 

e L (k;x,t) = PB R {-k\x,-t)P^ kedV L (5.40) 

as can be shown using the involutions ()3.22jl and ()5.39)1 . Substituting equations ()5.37)1 
and ()5.38jl into equation ()5.3()jl and using the results of subsection 13.41 one gets the 
following approximation of the matrix Or by a matrix Or 

0r = 0r[I + O(x- 1 -")] (5.41) 

uniform for k G &D R . Here 
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/ 



v R {k) 
where 



2 /r a k 1 



b R [x{k + 1)Y 



2A-„ 



2 /c o /c 1 



1-2 



• /j_\2 2k s 
1 \ 2ttJ 2fc»-fe-l 



!+5 \ 



/ 



(5.43) 



i7re ze 



2 -i(x-f) 



(5.44) 



lvre^e 1 ^ ^ 

It follows from the involution (|5.40|) that the matrix 9l is approximated by a matrix #l 



e L = h[i + o(x- 1 - e )] 

uniformly for k G &Dl, where 



9 L (k; x, t) = P9 R (-k; x, -t)p- x = I 



v L {k) 
x(k + 1) 



k e dV L 



and 



vi{k\ x, t) = —Pv R (—k; x, — t)P _1 
The explicit expression for vl reads 

/ • I n \2 2k s 



V L (k) - 

where 



2k s ~k+l 



i+5 



a L [x(l - fc)]" 



2/r.s 



2fe s -fc+l 



1--: 



2 A' s 



2/c s -fc+l 



i7re ae 



i§ -i(z+t) 



i7re ae 



iS„i(a;+t) 



(5.45) 



(5.46) 



(5.47) 



(5.48) 



(5.49) 



sin 2 f 

To calculate the conjugation matrix # s in the large x limit use the large x 
asymptotics of the integral in (|5.26|) 

dpe- T ^ G Q 



0{x 



-3/2\ 



(5.50) 



' s 27r p — k k — k s 

where Go is defined by equation ()2.19|) . Note that the leading term in (|5.5()j) is of the 
order of x~ x l 2 . Substituting equation ()5.50j) into equation ()5.26|) and using equation 
(JEH| one gets 

d s = d s [I + 0(x- 3 / 2 )] (5.51) 
uniformly for k G dT> s , where 



s (k)=I + 



v s {k) 



k G &D, 



and 



v s (k) 



-2G sin 2 UmY 




(5.52) 



(5.53) 



Zero temperature correlation functions for the impenetrable fermion gas 



30 



Consider the RH problem 

(a) S(k) is analytic in C \ (dV R U dV L U dV s ) 

(b) S+(k) = S-(k)9 a (k) kedV a a = R,L,s (5.54) 

(c) S(k) — > / as k — > oo 

with the jump matrices 0r, 6*^ and 9 S defined by (|5.42|) . ()5.46|) and (|5.52|) . respectively. 
In complete analogy with equation (J3.50)) the solution S(k) of the RH problem ()5.54j) 
approximates the solution S(k) of the RH problem ()5.29|) 

S{k) = S{k)[I + 0(x- 1 - 6 )] (5.55) 

uniformly outside a vicinity of the conjugation contours dT> R , &T>l and &D S . 

To solve the RH problem (|5.54j) represent S(k) as a product of two matrices 

S(k) = M{k)N{k) (5.56) 

where the matrix N(k) is the solution to the RH problem 

(a) N(k) is analytic in C \ (dV R U dV L ) 

(b) N+(k) = N-(k)6 R>L (k) ke&D R>L (5.57) 

(c) N(k) -> / as k -> oo 

with the jump matrices and ^ given by (|5.42j) and (|5.46|h respectively. Comparing 
(|5.54jl and ()5.57|) one sees that the matrix M(k) solves the following RH problem 

(a) M(k) is analytic in C \ dV s 

(b) M+(k) = M^NityOsityN-^k) k e dV s (5.58) 

(c) M(k) -»• / as fc -> oo 

where the jump matrix 6* s is defined by equation ()5.52|) . 

Let us solve the RH problem ()5.57|) first. Since the jump matrices (I5.42|) and 
()5.46|) have the same analytic properties in the complex fc-plane as the matrices ()3.45j) 
and (|3.48jl . the RH problem ()5.57|) can be solved similarly to the RH problem ([3.49)1 . 
Introduce the following notation for the analytic branches of N(k): 

( N R (k) keV R 
N(k) = I N L {k) keV L . (5.59) 

[ N O0 (k) keC\(V R UV L ) 

Using the same arguments as in deriving equation ()3.53|) from (J3.52)) one gets 

N 00 {k) = I+-^- + -^- keC (5.60) 
k — 1 k + 1 

where the matrices A R and Al do not depend on k. Dropping terms of order x~ 1 ~ e and 
higher, we find 

/ / 7r +'7 \ 
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and 

Al= 2_( -x l-itef^ ix L e^sin|e-(fc±i)"^^ 

" V ^ e " 13 csc ^ -*2 + £ (If) 2 t& , ' 



(5.62) 



Here 



2tv 



we"(2x)- l -i ( k] \— m-"(2x)- l +i ( k 2 

XR = pFPi I^FtJ " l = n(ft) 8 mi 1*FtJ (5 ' 63) 

and k is denned by equation ()3.64)1 . The matrices Nji(k) and N^k) are found from 

= ^(^[^(A;)]- 1 iV L (fc) = iVooC^fcCA:)]- 1 . (5.64) 



Next, solve the RH problem (J5.58jl . Use the same method as was exploited in 
solving the RH problems ()5.57|) and ()3.49|) . Denote the analytic branches of M(k) by 

The jump matrix NooOgN^ 1 in ()5.58l b) has a simple pole at k = k s . Therefore 

M OQ (k) = I + t— V" fceC (5.66) 

where the matrix A s does not depend on k. Choosing 

A s = N^Vsih^N^h)}- 1 (5.67) 
one solves (J5.58jl . This can be checked directly. 



5.6. Riemann- Hilbert problem \2.5°A) : results in the large x limit 

Like the t = case, the matrix K in the decomposition (J5.23)l does not contribute to 
the asymptotic solution of the RH problem (J5.22)) to the order we are interested in. 
Substituting (J5.55)l into the decomposition ()5.23)1 one gets the estimate 

U = SQ[I + 0(x-^ e )} (5.68) 

uniform outside an arbitrarily small vicinity of the contours &Dr^^ and S 13 s D T)^. 
Using ()5.15|) and ()5.2|) one gets the approximation to Y{k) up to the order of x~ 1 ~ e , 
uniform outside this vicinity. This completes the solution to the RH problem 1)2.52)1 in 
the large x limit. 



6. Time dependent correlation functions: space-like region 

In this section we derive the asymptotic expressions for the correlation functions Gh(x, t) 
and G e (x,t) in the space-like region |x| > 2\t\. The asymptotic expressions are derived 
assuming that x, t — > oo at a fixed value of the parameter k s defined in 1)5.1)1 . The 
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explicit expressions are written assuming x > and t > 0; all other cases follow from 
the relations fl2HJ) and (l2~2J) . 

All the calculations will be similar to those for the t = case. In subsect ion 16 . II we 
calculate the functions B ab and C ab f)2. from the large k expansion ()2.55|) of Y(k). 
Having these functions we calculate the Fredholm determinant det(J + V) in the large 
x limit, see equation (j6.8J) . The constant C(rj) entering ()6.8|) is calculated in subsection 
6.21 This constant is found to be the same as in the t = case. We give the answer for 
the correlation functions Gh{x,t) and G e (x, t) in subsection 16.31 



6.1. Large x asymptotics of B^, C ab and of the Fredholm determinant det(J + V) 

In this subsection we use the results of section on the asymptotic solution of the 
RH problem ()2.52|) to calculate the functions B ab , C a b and the Fredholm determinant 
det(J + V) in the large x limit. 

It follows from ijElEty . (J52D, jHHHD, flEZS) , fl£55J, and dHTEHD that in the 



(6.1) 



vicinity of k = oo the following uniform approximation is valid: 

Y = M 00 N 00 Q 00 V + 0(x~ 1 - e ). 



Recall that the matrix is given by equation 



the matrix iVoo by equation 



()5.60|) . the matrix by ()3.21|) . the matrix V by ()5.3|) . The approximation (jSH} being 
uniform, the functions B a b and ()2.16j) can be calculated in the large x limit from 
the large k expansion ()2.55|) of equation (|6.1|) . In particular, 



= 2ie i 'e^^ esc- 



"~2tT 



XL 



2ir 



[1 + 0{x~ 



(6.2) 



The function entering (J2.22j) has the form 



2ie : *e 2 sin - 



1 
2 



fee 



7T — 



XL 



fc, - 1 



2ir 



— 2sm - 



k s 

2 V /c .e 



Go HI + 0( 



x 



(6.3) 



The expressions for the other potentials are quite bulky, so we write down the differential 
equations (j2.41|) and (|2.42j) in the large x limit omitting intermediate calculations 



a x lndet(/ + V) ~- 



it] 

7T 



ki 



Tf 



k 2 -l 2tt 2 x 



+ 2in' 1 (>c 2 R - 



x7 



1)k + 2 % + 1)k-2 



and 



(A* - l)K 

a t lndet(/ + t>) 



2k $ 



1)k 



a L + (x- 1 -*) 



T} 



k 2 -l 2n 2 x 



(6.4) 



(6.5) 
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where 

/ /■ _ 1 \ 



a R = 2ie 1 '(l - e-^^oK -1 ( ) (6.6) 



3(7T-Tj) 
fc, - 1 \ 2 - 



a L = -2ie"(l - e-^^GoK- 1 J (6-7) 

and k is given by equation (|3.64l) . 

Integrating equations ()6.4|) and ()6.5|) asymptotically one gets 

In det(J + V) = ^x - (J^y In [x 2 (l - A;; 2 )] + C(rj) + O (x' e ) . (6.8) 

Note that x 2 {l — k~ 2 ) = xrXl- We calculate the integration constant C(rj) in subsection 
16.21 It will be the same as for the t = case. 

6.2. Calculation ofC(r]) 

In this subsection we calculate the constant C{ri) in (|6.8|) using the differential equation 
(I2.37|) . The calculation will resemble that of subsection I4.21 

Consider the first term in the right hand side of equation f)2.37j) . Split the integral 
in equation (|2.38j) in two parts and define 

Jr = r^— I dk[f(k)fa 2 d k f(k) (6.9) 
1 — e Jo 

Jl = J dk[f(k)fa 2 d k f(k). (6.10) 

Consider Jr first. The large x asymptotics of the function f(k) in equation ()6.9j) is 
calculated as follows. Use equation ()2.53|) 

f(k) = Y + (k)e(k) k G S 2 . (6.11) 

Substituting (Q and (JHHSJ) into (jHHID and using (j2,llj) and (j2.25j) one gets 

f(k) = Z + (k)e (k) keY, 2 (6.12) 

where 

the matrix Z(&;) solves the RH problem ()5.10|) and the functions e„(/c) and r(k) are 
defined by equation (|2.11|) . It follows from f|5. 15|) and (|5.16|) that 

Z+(k) = U + (k)m(k) k G S 2 (6.14) 

where £/(&) solves the RH problem f)5.22j) . Substituting equation ()6.14|) into ()6.1H) and 
performing the analytic continuation of the formula (|5.68|) in complete analogy with 
the analytic continuation of the formula ()3.67|) discussed in subsection 14.21 and shown 
in figure El we obtain the approximation to f(k) 

f = M^NrQrII + C^aT 1 - 6 )] (6.15) 
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uniform for k G [0,1]. The matrices and N R are defined in subsection 15.51 the 
function g is given by the formula 

g R (k) 

Substituting into (jfj.lfjjl the expression for Q R (k), defined by equations (|5.37|) . ()5.33j) 
and ()3.31|) . and using the identity ()4.13|) one gets 

R(u ^z T{k)/2 - T{1) 1F1 [& + 1, i; ten&R -i)]fk-i\* 

g+[k;x,t) - 1 1 



Mk) = I y t) u [ ) = [Q*(A:)] W*0eoW fc € [0, 1]. (6.16) 



9-(k;x,t) 



R 



(6.17) 

__ 1 \ 2. 

rfe)sm2 [**(* + 1)]"* VAh"^TJ 

where A# = A j r(/c) is defined by equation (|5.32j) . 

Using the uniform estimate ()6.15|) and equations ([6.17)1 . ()5.60|) - ([5.67p one can derive 
the following estimate for ()6.9|) 

Jr = t^— ! dkigni^fa^gnim + Oix- 1 -')}. (6.18) 
1 - e vr] Jo 

Performing the integral in ()6.18|) asymptotically one gets 

ix R T) 1 

n mx R -\ — i 

2tt 2tt 2 K 2 

where Ci{rf) is defined by equation (|4.17j) . 

The large x asymptotics of Jl is calculated similarly, yielding 

Jl = ^ " 2^ lnXL + ^ l(r?) + ° {x ' % (6 ' 20) 
We have calculated the large x asymptotics of the first term in the right hand side 
of equation (|2.37p . By analysis of equation ()2.40|) one can show that the second term 
in the right hand side of equation (|2.37j) is of the order of x~ e . Therefore, the large x 
asymptotics of (J2.37)) is 

d v In det(f + V) = J R + J L + 0(x~ e ) 

= --^L \n(x R x L ) + CM + <D(x-*). (6.21) 

Integrating equation ()6.21|) over 77, taking into account the initial condition (|2.43|) 
and comparing the resulting expression with ()6.8|) one finds that C{rf) is given by the 
same expression (|4.19jl as in the t — case. 



J R = ^ ~ lnx * + H C M + £>(x~ e ) (6-19) 



6.3. Time dependent correlation functions in the space-like region: the results 

In this subsection we write down the asymptotic expressions for the correlation functions 
Gh{x,t) and G e (x,t) in the space-like region \x\ > 2\t\. The asymptotic expressions are 
given assuming that x, t — > 00 at a fixed value of the parameter k s defined in (15. lj) . The 
explicit expressions are written assuming x > and t > 0; all other cases follow from 
the relations fl2IEJ) and 
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The correlation function Gh{x,t) is given by the formula (|2.32|) . The asymptotics 

of B is given by ()6.2}) ; the asymptotics of det(J + V) by ()6.8|) with the constant C{rf) 

given by (j4.19|) . Recall that z = exp(i?7), therefore the point z = 1/2 corresponds to 
77 = iln2. The leading contribution to the asymptotics of Gh{x,t) comes from the first 
term in the right hand side of equation ()2.32|) . yielding 

G h (x,t) 



cxp 



2i (x + 2t) 2A (x-2t) 2A 

^ g— a;ln2/7r 

_ 2i(x-2t) 2A (x + 2t) 2A 



iln2 n , 

i(^x — x J — m(x + It) 



Ti 



exp 



-IX 



x . 



iln2, . 

H m(x 

7T 



2t) 



(6.22) 



where the normalization constant S is defined by equation ()4.23|) and the phase shift xo 
by equation ()4.22|) . The anomalous exponents are given by 

2 1 /, o\ 2 



A 



1 



In 2 

7T 



A 



In 2 



7T 



(6.23) 



The relative correction to ()6.22|) is of the order of x~ l . 

Show that the second term in the right hand side of equation (J2.32j) does not 
contribute to the asymptotics ()6.22|) . By virtue of ()6.8|) . the determinant det(7 + V") is 
of the order of exp(xlnr/7r) on the integration contour 7, where 7 is the circle \z\ = r, 
figure Therefore, for r < 1/2 and for sufficiently large x the contribution of the second 
term is negligible. 

The correlation function G e (x,t) is calculated similarly, yielding 



G e (x,t) = -G* h (x,t) +e it G (x,t) 



e -xln2/7r e C(iln2) 



x-2t 



iln2/7T 



(6.24) 



(x 2 -4t 2 ) 2A \x + 2t / 
where Gh(x,t) is given by equation ()6.22|) . the function Go(x,t) is defined by equation 
(g3Sl) and A defined by fl6~2H . 

Compare the result ()6.24|) with the general relation ()2.3|) . One can see that 
the last term in the right hand side of equation 1)6.24)1 corresponds to the averaged 
anticommutator of fermion fields in equation ()2.3)) . It contains a rapidly oscillating 
factor exp (it). This means that although formally present in the long distance expansion 
of the correlation function G e this term is irrelevant at energies smaller than the chemical 
potential. 



7. Riemann-Hilbert problem at t ^ 0. Time-like region. 

In this section we find the asymptotic solution of the Riemann-Hilbert problem (J2.52)) 
for |x| < 2\t\. We assume that x > and t > 0. We will consider the large t asymptotics 
at a fixed value of the parameter k s defined in ()5.1|) . 

In subsection 17. II we reformulate the RH problem ()2.52)) so as to make it similar to 
the RH problem ()3.3|) . In subsection 7.2 we repeat the construction made for the t = 
case in subsections 13.21 and 13.31 In subsection 7.3 we give an explicit solution to the 
jump relation (7.7.b) in a vicinity of the points k = ±1 and k = k s . In subsection 7.4 
we perform the matching procedure analogous to the one of subsection 13.51 
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Figure 11. The conjugation contour Si U S 2 U S3 U S 4 for the RH problem (|7.7|l . 
The contour Si consists of two semi-infinite lines parallel to the imaginary axis. 



7.1. Reformulation of the Riemann-Hilbert problem 

In this subsection we reformulate the RH problem (J2.52)) so as to make use of the results 
of subsections 15.31 and 15.41 

Consider a contour shown in figure ^3 This contour divides the complex /c-plane 
into four regions, denoted as S\ through S\y. Represent the solution Y{k) of the RH 
problem (p2~5~2|) product of two matrices 

Y{k) = Z{k)V{k) (7.1) 

where 

v(k) _( 1 \( 1 -2sin 2 f[£(A;) -ie~^o 2 (k)] \ 

V[k) - [2ie^ gi (k) 1 ){o 1 J' {7 - 2) 

The piecewise constant functions Q\ and Q2 are defined by 

f e ir > ke Sa 

Qi{k) = I -e' ir > k e S w (7.3) 
I otherwise 

and 

{1 k G S m \mk > 

-1 k G Si Imk < . (7.4) 
otherwise 

The contour Si is chosen so that the real part of r(k) is positive for all k G Si except 
k — ±1 and 

Rer(A;) — » +00 as k — > 00 G Si (7.5) 

Note that the condition (|7.5|) implies 

V(k) — > I as fc -> 00. (7.6) 

Using equations ()2.52|) . (|7.1|) and ()7.2|) . and the property ()7.6|) one gets the following 
RH problem for Z{k) 

(a) Z(k) is analytic in C \ S 

(b) Z + (k) = Z_(k)ni(k) k G Si i= 1,2,3,4 (7.7) 

(c) Z(k) —>■ I as k — > 00 
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Figure 12. The contour used in the factorization of the RH problem (|7.7I) . The discs 
T^l and 2? s are shaded in gray. 

where the jump matrices /jli, /U 2 and /i 3 are given by equations ()5.17|1 . (|5.2(Jj) and (|5.18|) . 
respectively. The matrix /i 4 is defined as 

^ = ( 2ie-«e i " 1 ) • (7 ' 8) 
7.2. Factorization of the Riemann-Hilbert problem \7. 7p 

In this subsection we sketch the scheme of the analysis of the RH problem f)7.7j) in the 
large t limit. This scheme is similar to the one given in subsection 15.31 

The matrices fi±, fi^ and /14 converge to the identity matrix, as t — * 00, everywhere 
on their conjugation contours except at the points k = ±1 and k = k s . Introduce three 
disks V L , V R and V s as shown in figure IT21 Denote the domain C \ (V L U V R UV S ) as 
Doq. Represent the solution Z(k) of the RH problem ()7.7|) as a product of three matrices 

Z = KSQ. (7.9) 

Here the matrix Q(k) is defined by 

Q(k) = Q a (k) keV a a = L,R,s, 00. (7.10) 

The matrix Qoo is given by equation ()3.18|) . The matrices Qr and Ql are defined in Vr 
and Vl, respectively, and satisfy there the jump relations 

[Q R (k)} + = [Q R (k)}_ f i i (k) ke^nVn 2 = 1,2,4 (7.11) 

and 

[Q L (k)] + = [Q L (k)]_ f i l (k) ke^nv L i = 1,2,3. (7.12) 

The matrix Q s {k) is defined in T> s and satisfies 

[Q s {k)} + = [Qs(k)]_tn(k) keE { nv s i = 2,3,4. (7.13) 

We give solutions to the jump relations ()7.11|) . ()7.12|) and ()7.13|) in subsection 17.31 

The matrix S(k) in the decomposition ()7.9|) is the solution to the RH problem 
(|5.29|) with the conjugation matrices 9r^ jS defined by equation (|5.3U|) . The matrices 
Qr, Ql and Q s in ()5.30|) solve the jump relations (j7.11|) . f)7. 12j) and ()7.13|) . respectively. 
Like the spacelike case, the matrix K can be replaced by the identity matrix. 
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7.3. Exact solutions to the jump relations J7. 1 J7. 1H\} and \7.1ty 

In this subsection we solve explicitly the jump relations (|7.11jl . (|7.12|) and (|7.13jl . We 
use the results of subsection 15.41 

The jump relation ()7.12|) coincides with the jump relation ()5.28|) . Therefore one 
can use the solution 1)5.39)1 to this relation. To obtain a solution to the jump relation 
(j7.11j) consider the conjugation matrices fa, fa and fa defined on the contours S 1; S 2 
and £4, shown in figure El These matrices satisfy the following involutions 



fa(k;x,t,r}) =p- 1 fa(k 
l^2(k;r]) =P{ l fa(k 
fa(k;x,t,r]) =P~ 1 fa(k 



-x, -t, -i])P 1 

-v)Pi (7-14) 
-x, -t, -rj)Pi 

where P\ is a ^-independent matrix defined by 

pi -ie-W sin 2 \ 

Pl -^ie^ 2 csc| )■ (7 ' 15) 

Denote the solution ()5.37|) to the jump relation ()5.27|) by Q%(k; x, t, rj). Comparing the 
conjugation contours in the disks T> R shown in figures El and ^] and using the involutions 
(17.14)) a solution Q R to the jump relation ()7.11j) can be found from 

Q R (k; x, t, 77) = PiQnik; -x, -t, -^P^ 1 (7.16) 

Finally, write down a solution to the jump relation ()7.13j) : 

Qs(k)= I I " }Qoo{k) (7.17) 



q(k) 



where 



and the integral runs along the contour S s which consists of the part of the contour S 3 
lying outside V L and the part of the contour S 4 lying outside V R , see figure IT21 

7-4- Solution to the Riemann-Hilbert problem in the large t limit 

In this subsection we solve the RH problem ()5.29j) in the large t limit. The procedure 
will resemble that of subsection 15.51 

The conjugation matrix 9l is approximated by a matrix 9l 

e L = ~e L + o (r x - e ) (7.19) 

uniformly for k G OVl, where 9l is defined by equations (|5.4fij) and (|5.48j) . The 
conjugation matrix 6 R is approximated by a matrix 6 R 

9r = R + O (r 1 - 6 ) (7.20) 

uniformly for k G &D R , where the matrix 9 R is given by equation ()5.42)1 and v R is defined 

as 

v R (k; x, t, 77) = -Pivi(k; -x, -t, -^P^ (7.21) 
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Here we use the notation v s ^(k; x, t, rj) for the matrix given by equation (|5.43|) . The 
matrix Pi is defined by equation (j7.15|) . The conjugation matrix 9 S is approximated by 
a matrix 9 S 

9 s = 9 s [I + 0(t^ 2 )] (7.22) 

uniformly for k G &D S . The matrix 9 S is given by equation ()5.52|) where v s (k) is 
calculated using (fOHjl . (|7T7j) and (fTTHjl 







i-fc, 
i+fc 6 







G* 



(7.23) 



The asymptotic estimate for the matrix S is given by equation ()5.55|) where the 
matrix 5* is the solution of the RH problem ()5.54j) with the matrices 9r^,s defined in 
the paragraph above. Like the spacelike case, solve the RH problem (|5.54|) using the 
decomposition f)5.56j) . The matrix entering the definition ()5.59|) is given by equation 
(I5.6(J|) with matrices A R and A L given by 



.4 



R 



and 



.4, 



Here 



2 

k 



/ 



l-k. 



TV — 7] 



-x R sin 2 e - i4 



^2 _ j_ ( V\ 2 1 
■R 4t V27T/ 1-* 



+■7 \ 
2tt 1 



-l-e> 



(7.24) 



\ -^cscfe 1 * 



7T + ?J 

2tt 



, 1 (j_\ 2 _J_ 

~ it \2tt) 1+k 



+ 0(f 



(7.25) 



XL 



me lx {Aty 



r2 {-£) e ' 2 ' sin 2 

ivre- ia; (4t)- 1+ - 



1 



2 V 1 - k 2 
1 



2tt 



(7.26) 



and k is defined by equation ()3.64|) . The matrix M OQ (k) in the decomposition ()5.56|) is 
given by equation (|5.66|) with A s defined by equation (|5.67|) . 



8. Time dependent correlation function: time-like region, k s ^ 

In this section we derive the asymptotics of the correlation functions G e (x, t) and Gh(x, t) 
in the time-like region \x\ < 2\t\. The asymptotic expressions will be given assuming 
that x, t — > oo at a fixed value of the parameter k s defined in equation (|5.1|) . It will also 
be assumed that x > and t > 0; all other cases follow from the relations (j2.1|) and 

(JZ21). 

The calculations in this section are similar to those given in sections 0] and |H1 
We will therefore omit most technical details already explained in those sections and 
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focus on the particulars of the time-like case. In subsection 18.11 we give the the large 

t asymptotics of B , b ++ and det(7 + V). In subsection 18.21 we perform the contour 

integral entering representations ()2.21|) and ()2.22j) and obtain the asymptotic expressions 
for the correlation functions G e (x,t) and Gh{x,t). 



8.1. Large t asymptotics of B^, C a b and of the Fredholm determinant det(J + V) 

The asymptotic expressions for the functions B and b ++ are calculated similarly to 

subsection 16.11 yielding 



B_ 



1 - k. 



o 



2e i4 esc - 



1 + k. 



TV — 11 

~~2tT 



xl + 



1 - k s 
l + k s 



rr-r-y 
2ir 



[1 + 0(t~ 



and 



2e 



sm ■ 



l-h 



7T — I] 



XL 



l + k s 



rv + y 
2tt 



[i + o(r 



•i) 



1.2) 



where xr and are given by equaiton ()7.26)1 . The large t asymptotics of the Fredholm 
determinant det(J + V) is 



mdet(J + y) 



7T 



-X — 



2ll 



In [At 2 (l - kf)] + C( V ) + O (r e ) . (8.3) 



The constant C(r]) can, in principle, be calculated in the same manner as the 
constant C(ri) for the space-like region. There is, however, a difference which did not 
allow us to obtain a closed analytical expression for C(rj) so far. Unlike the space-like 
case, the second term in the right hand side of equation ()2.37|) is not small as a function of 
the asymptotic parameter t. Therefore, the calculation of C (if) requires the asymptotic 
analysis of the expression ()2.40|) . which is straightforward but cumbersome. This is in 
contrast with the high temperature case, where the calculation of such a constant was 
done using an additional differential equation involving the functions B a b and C a t only 
[TT] . Such an equation does not exist in the zero temperature case, considered here. 



8.2. Time dependent correlation functions: time-like region, k s ^ 0. The results 

In this subsection we write down the asymptotic expressions for the correlation functions 
Gh{x,t) and G e (x, t) in the time-like region \x\ < 2\t\. The asymptotic expressions are 
given assuming that x, t —>■ oo at a fixed value of the parameter k s defined in (I5.1J1 . The 
explicit expressions are written assuming x > and t > 0; all other cases follow from 
the relations fl2HJ) and (jZfl . 

Consider the representations ()2.32|) and ()2.33|) . Applying the same arguments as 
given in subsection (jfi.3)l one can show that the contributions from the integrals along 
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the contour 7 in ([2.32)1 and ([2.33)1 can be neglected in the large t limit. Substituting 
the results flEU, (Q and (Q in equations (l2~32j) and (l2~331) we find 



G h (x,t) 



and 

G e (x,t) 



-x In 2/tt 



2i (2t + x) 2A {2t - x) 2A 



exp 



iln2 



1{X — Xo) 



71 



ln(2t + x) 



-x In 2/-7T 



+ 



4i (2t-x) 2A (2t + x) 2A 

i-f 



exp 



iln2, . . 
— \[x — x ) H m(2t — x) 



7T 



e» -xln2A e C(iln2) / 2 t-a; 



-i In 2/tt 



-x In 2/-7T 



i (2t-x) 2A (2t + x) 2A 

S g— a;ln2/7r 

" 2i(2t + x) 2A (2t-a;) 2A 



exp 



i In 2 . 
i(x — xq) m(2t — x) 



7T 



exp 



i In 2 , . 
-i(x — xo) H — m(2t + x) 



7T 



(8.4) 



i.5) 



where the dimensions A and A are given by equation (|f).2H|l . the phase shift xo by 
equation ([4.22)1 . The normalization constant H can be obtained from the constant S 
given in equation ()4.23|) by replacing the constant C{rf) with C(r)). 



9. Time-dependent correlation functions: time-like region, x = 

In this section we consider the special case of x = 0. The large t asymptotics of the 
correlation functions G e (0,t) and Gh{0,t) are given assuming that t > 0; the case 
of negative t follows from the relation (12.2(1 . In the derivation of the asymptotical 
expressions we the use the results of section |HJ 

For x = it is convenient to calculate the large t asymptotics of the integrals 
(2.6) and (2.7) without the deformation of the integration contour leading to (2.32) and 
(2.33). Substituting asymptotic formulae (8.2) and (8.3) in equation (2.7) we arrive at 
the expression 



G e (0,t) 



drj 



F(rj) i4" 



+ 



l-cosr? r 2 (^)e- 



exp 



-1 



1 

TV 



2tt 



hit 



dr]- 



exp 



V 



-1- J--2 p- 

7T V27T 



V 



lnt \ . (9.1) 



1-cost] r 2 (-^)e 1 i 

For large t the integrals on the right hand side of (9.1) can be evaluated in the saddle 
point approximation with lnt playing the role of the asymptotic parameter. The relative 
correction to the integrands in equation (9.1) is of the order of t~ e . At the saddle points, 
7] = ±7T, one has e = leading to an unknown constant factor. This factor can, 
in principle, be calculated using the approach presented in this paper. The resulting 
expression for G e reads 



G e (0,t) 



const 
Vt hit 



1 + 



(9.2) 
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The integral (2.6) is evaluated similarly, yielding 

G h{O , t) = ^L + A-"G-(0,t) (9.3) 
Vtlnt Vint 

with the relative correction of the order of 1/ Int. 
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